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(57) ABSTRACT 

Disclosed herein are systems, computer-implemented meth- 
ods, and tangible computer-readable storage media for wide 
field imaging interferometry. The method includes for each 
point in a two dimensional detector array over a field of view 
of an image: gathering a first interferogram from a first detec- 
tor and a second interferogram from a second detector, modu- 
lating a path-length for a signal from an image associated with 
the first interferogram in the first detector, overlaying first 
data from the modulated first detector and second data from 
the second detector, and tracking the modulating at every 
point in a two dimensional detector array comprising the first 
detector and the second detector over a field of view for the 
image. The method then generates a widefield data cube 
based on the overlaid first data and second data for each point. 
The method can generate an image from the widefield data 
cube. 

19 Claims, 14 Drawing Sheets 
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SYSTEMS, COMPUTER-IMPLEMENTED 
METHODS, AND TANGIBLE 
COMPUTER-READABLE STORAGE MEDIA 
FOR WIDE-FIELD INTERFEROMETRY 

5 

ORIGIN 

The disclosure described herein was developed by employ- 
ees of the United States Government and may be manufac- 
tured and used by or for the United States Government for 10 
governmental purposes without the payment of any royalties 
thereon or therefore. 

BACKGROUND 

15 

1 . Technical Field 

The present invention relates to interferometry and more 
specifically to wide field interferometry. 

2. Introduction 

Michel son imaging interferometry only allows for very 20 
narrow field-of-view imaging and spectroscopy. It is gener- 
ally accepted that the field-of-view is limited to the angular 
extent of the core of the primary beam of ±X/D, where X is the 
wavelength and D is the diameter of a single aperture. How- 
ever this limitation is not a fundamental physical limitation 25 
but arises due to the Van Cittert-Zemike theorem and single 
pixel or small format detectors. This interferometry approach 
is used in observing and measuring in space to construct 
images from multiple smaller, inexpensive telescopes, which 
rival images generated from much larger, more expensive 30 
telescopes which are difficult to construct and maintain. How- 
ever, this traditional interferometry approach generates very 
narrow field images. Accordingly, what is needed in the art is 
improved wide field imaging. 

35 

SUMMARY 

Additional features and advantages of the disclosure will 
be set forth in the description which follows, and in part will 
be obvious from the description, or can be learned by practice 40 
of the herein disclosed principles. The features and advan- 
tages of the disclosure can be realized and obtained by means 
of the instruments and combinations particularly pointed out 
in the appended claims. These and other features of the dis- 
closure will become more fully apparent from the following 45 
description and appended claims, or can be learned by the 
practice of the principles set forth herein. 

Disclosed are systems, computer-implemented methods, 
and tangible computer-readable storage media for wide field 
imaging interferometry. The method includes for each point 50 
in a two dimensional detector array over a field of view of an 
image: gathering a first interferogram from a first detector and 
a second interferogram from a second detector, modulating a 
path-length for a signal from an image associated with the 
first interferogram in the first detector, overlaying first data 55 
from the modulated first detector and second data from the 
second detector, and tracking the modulating at every point in 
a two dimensional detector array comprising the first detector 
and the second detector over a field of view for the image. The 
method further generates a widefield data cube based on the 60 
overlaid first data and second data for each point. The method 
can generate an image from the widefield data cube. The first 
and second detectors can be selected from a group of three or 
more detectors. Either the first and/or second detectors can be 
different physical detectors for each point. The method can 65 
further gather more than one first and second interferograms 
at different distances between the first and second detectors 
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for each point in the two dimensional detector array. The 
method can gather more than one first and second interfero- 
grams at different locations for the first and second detectors 
for each point in the two dimensional detector array. The 
interferograms can represent visible light, non-visible light, 
sound, or other waves. In one embodiment, the detectors are 
telescopes. 

BRIEF DESCRIPTION OF THE DRAWINGS 

In order to describe the manner in which the above-recited 
and other advantages and features of the disclosure can be 
obtained, a more particular description of the principles 
briefly described above will be rendered by reference to spe- 
cific embodiments thereof which are illustrated in the 
appended drawings. Understanding that these drawings 
depict only exemplary embodiments of the disclosure and are 
not therefore to be considered to be limiting of its scope, the 
principles herein are described and explained with additional 
specificity and detail through the use of the accompanying 
drawings in which: 

FIG. 1 illustrates an example system embodiment; 

FIG. 2 illustrates an example simplified schematic of an 
interferometer; 

FIG. 3 illustrates an example wide-field imaging interfer- 
ometry testbed (WIIT); 

FIG. 4 illustrates an example data format for WIIT data; 

FIG. 5 illustrates example WIIT interferograms and vis- 
ibilities; 

FIG. 6 illustrates an example WIIT fractal bar target and 
images; 

FIG. 7 illustrates an example image of WIIT one-dimen- 
sional ultraviolet plane coverage; 

FIG. 8 illustrates an exemplary WIIT dirty beam; 

FIG. 9 illustrates example WIIT spectra for varying base- 
lines; 

FIG. 10 illustrates sample WIIT dirty image data results 
and maximum entropy deconvolutions; 

FIG. 11 illustrates an example inverse Fourier transform 
over B with 0 as a conjugate variable; 

FIG. 12 illustrates an exemplary method embodiment; 

FIG. 13 illustrates an example configuration of two detec- 
tors on a boom; 

FIG. 14 illustrates an example configuration of two detec- 
tors on a rotatable boom; and 

FIG. 15 illustrates an example grid configuration of more 
than two detectors. 

DETAILED DESCRIPTION 

Wide-field imaging interferometry testbed (WIIT) can 
demonstrate imaging and spectroscopy over fields-of-view 
larger than the typically narrow primary beam footprint of a 
conventional Michel son interferometer. At the heart of this 
technique is the “double-Fourier” approach in which the aper- 
tures and a delay line are both moved to collect interfero- 
grams over a two dimensional wide field detector grid simul- 
taneously, one interferogram per detector pixel. This 
aggregate set of interferograms, as a function of baseline and 
delay line, is algorithmically processed to construct a hyper- 
spectral image cube. Algorithms and approaches to construct 
the image cube are disclosed herein. 

Classical Michelson imaging interferometry allows for 
high spatial and spectral resolution imaging over a small field 
of view (FOV). The double-Fourier Michelson technique 
expands the FOV, allowing for wide field imaging and spec- 
troscopy of each pixel of the high resolution, wide-field 
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image. In double Fourier wide field imaging interferometry, 
the fields are combined in the pupil plane but the data is 
collected in a wide format focal plane. A baseline algorithm 
converts the collected data to a hyperspectral data cube con- 
sisting of two dimensional high resolution wide format 
images, one image per wavelength (or frequency). A detailed 
mathematical analysis, set forth below, includes polarization 
and other effects. A physical optics model can estimate the 
dirty hyperspectral image data cube and ultimately deconvo- 
lute the estimated image cubes with missing or attenuated 
spatial frequencies. The treatment is via scalar field theory 
with no polarization effects. One advantage of this approach 
is that the construction of the dataset and the wide-field imag- 
ing algorithms can be concisely understood with the equa- 
tions and derivations set forth below. 

As a preliminary matter, various embodiments of the dis- 
closure are discussed in detail below. While specific imple- 
mentations are discussed, it should be understood that this is 
done for illustration purposes only. A person skilled in the 
relevant art will recognize that other components and con- 
figurations may be used without parting from the spirit and 
scope of the disclosure. 

With reference to FIG. 1, an exemplary system 100 
includes a general-purpose computing device 100, including 
a processing unit (CPU or processor) 120 and a system bus 
110 that couples various system components including the 
system memory 130 such as read only memory (ROM) 140 
and random access memory (RAM) 150 to the processor 120. 
These and other modules can be configured to control the 
processor 120 to perform various actions. Other system 
memory 130 may be available for use as well. It can be 
appreciated that the disclosure may operate on a computing 
device 100 with more than one processor 120 or on a group or 
cluster of computing devices networked together to provide 
greater processing capability. The processor 120 can include 
any general purpose processor and a hardware module or 
software module, such as module 1 162, module 2 164, and 
module 3 166 stored in storage device 160, configured to 
control the processor 120 as well as a special-purpose pro- 
cessor where software instructions are incorporated into the 
actual processor design. The processor 120 may essentially 
be a completely self-contained computing system, containing 
multiple cores or processors, a bus, memory controller, cache, 
etc . A multi-core processor may be symmetric or asymmetric . 

The system bus 110 may be any of several types of bus 
structures including a memory bus or memory controller, a 
peripheral bus, and a local bus using any of a variety of bus 
architectures. A basic input/output (BIOS) stored in ROM 
140 or the like, may provide the basic routine that helps to 
transfer information between elements within the computing 
device 100, such as during start-up. The computing device 
100 further includes storage devices 160 such as a hard disk 
drive, a magnetic disk drive, an optical disk drive, tape drive 
or the like. The storage device 160 can include software 
modules 162, 164, 166 for controlling the processor 120. 
Other hardware or software modules are contemplated. The 
storage device 160 is connected to the system bus 110 by a 
drive interface. The drives and the associated computer read- 
able storage media provide nonvolatile storage of computer 
readable instructions, data structures, program modules and 
other data for the computing device 100. In one aspect, a 
hardware module that performs a particular function includes 
the software component stored in a tangible and/or intangible 
computer-readable medium in connection with the necessary 
hardware components, such as the processor 120, bus 110, 
display 170, and so forth, to carry out the function. The basic 
components are known to those of skill in the art and appro- 
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priate variations are contemplated depending on the type of 
device, such as whether the device 100 is a small, handheld 
computing device, a desktop computer, or a computer server. 

Although the exemplary embodiment described herein 
5 employs the hard disk 160, it should be appreciated by those 
skilled in the art that other types of computer readable media 
which can store data that are accessible by a computer, such as 
magnetic cassettes, flash memory cards, digital versatile 
disks, cartridges, random access memories (RAMs) 150, read 
to only memory (ROM) 140, a cable or wireless signal contain- 
ing a bit stream and the like, may also be used in the exem- 
plary operating environment. Tangible computer-readable 
storage media expressly exclude media such as energy, carrier 
signals, electromagnetic waves, and signals per se. 

15 To enable user interaction with the computing device 100, 
an input device 190 represents any number of input mecha- 
nisms, such as a microphone for speech, a touch-sensitive 
screen for gesture or graphical input, keyboard, mouse, 
motion input, speech and so forth. The input device 190 may 
20 be used by the presenter to indicate the beginning of a speech 
search query. An output device 170 can also be one or more of 
a number of output mechanisms known to those of skill in the 
art. In some instances, multimodal systems enable a user to 
provide multiple types of input to communicate with the 
25 computing device 100. The communications interface 180 
generally governs and manages the user input and system 
output. There is no restriction on operating on any particular 
hardware arrangement and therefore the basic features here 
may easily be substituted for improved hardware or firmware 
30 arrangements as they are developed. 

For clarity of explanation, the illustrative system embodi- 
ment is presented as comprising individual functional blocks 
including functional blocks labeled as a “processor” or pro- 
cessor 120. The functions these blocks represent may be 
35 provided through the use of either shared or dedicated hard- 
ware, including, but not limited to, hardware capable of 
executing software and hardware, such as a processor 120, 
that is purpose-built to operate as an equivalent to software 
executing on a general purpose processor. For example the 
40 functions of one or more processors presented in FIG. 1 may 
be provided by a single shared processor or multiple proces- 
sors. (Use of the term “processor” should not be construed to 
refer exclusively to hardware capable of executing software.) 
Illustrative embodiments may comprise microprocessor and/ 
45 or digital signal processor (DSP) hardware, read-only 
memory (ROM) 140 for storing software performing the 
operations discussed below, and random access memory 
(RAM) 150 for storing results. Very large scale integration 
(VLSI) hardware embodiments, as well as custom VLSI cir- 
50 cuitry in combination with a general purpose DSP circuit, 
may also be provided. 

The logical operations of the various embodiments are 
implemented as: (1) a sequence of computer implemented 
steps, operations, or procedures running on a programmable 
55 circuit within a general use computer, (2) a sequence of com- 
puter implemented steps, operations, or procedures running 
on a specific -use programmable circuit; and/or (3) intercon- 
nected machine modules or program engines within the pro- 
grammable circuits. The system 100 shown in FIG. 1 can 
60 practice all or part of the recited methods, can be a part of the 
recited systems, and/or can operate according to instructions 
in the recited tangible computer-readable storage media. 
Generally speaking, such logical operations can be imple- 
mented as modules configured to control the processor 120 to 
65 perform particular functions according to the programming 
of the module. For example, FIG. 1 illustrates three modules 
Modi 162, Mod2 164 and Mod3 166 which are modules 
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configured to control the processor 120. These modules may 
be stored on the storage device 160 and loaded into RAM 150 
or memory 130 at runtime or may be stored as would be 
known in the art in other computer-readable memory loca- 
tions. 

Having disclosed some basic system components, the dis- 
closure now turns to the principles of wide-field Michelson 
imaging and specifically the mathematical model of wide- 
held Michelson imaging. The method embodiments set forth 
herein are discussed in terms of systems, such as the system 
100 shown in FIG. 1, configured to practice the method. An 
interferometer operates based on the physical optics model of 
propagation through the interferometer to generate a focal 
plane image. FIG. 2 illustrates a simplified example sche- 
matic of an interferometer 200 for reference. Equation 1, 
shown below, describes the scalar electric held, at the plane of 
the interferometer’s entrance pupil 202, due to a given emis- 
sion point 204 on a distant source at a single optical fre- 
quency: 


E 0 (a, fa, V) = \E 0 (U, fa, V )\ e ma^v) e i-—{ax + py) e i—z e -i^vt 


(2) 


This represents the equation of a wave arriving from a sky 
angle of (a,p) and at frequency v. The total held at point 
(x,y,z) and time t, from all source frequencies and held angles, 
is given by the integral over all sky angles in the held of view 
and frequencies in the pas sb and as: 


E 0 (x, y, z; t ) = 


(3) 
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fff I Eo(a, P v)| 


which takes the form of a 3D Fourier transform of IE 0 (a,(3; 
v)| e ^( a, ^ J v) and where 


20 


z 

c 


t 


E 0 (a, fa v) = \E 0 (a, fa, v)|e W(a ^ v >e^V“ ,rw 

where (a,|3) are the sky angles, i.e. angular offsets from the 
line of sight pointing of the interferometer, and thus (a, (3)= 
(0,0) dehnes the center of the held of view of the interferom- 
eter, i.e. its pointing. The optical frequency is given by v=c/X 
where c is the vacuum speed of light and X the vacuum 

wavelength, r =(x,y,z) are coordinates in the entrance pupil 
plane 202, where (x,y) are coordinates along the plane of the 
aperture and z is the distance from the source 204 to the 
aperture. IE 0 (a,®;v)l is a real valued function representing 
the held amplitude of a source point and a given frequency. 
0(a,|3;v) represents a uniformly distributed random phase on 
the interval of [-jt,jt] due to each source point at each fre- 
quency having a phase which is uncorrelated with all others 
over sky angles and frequencies. Each distant source point 
204 emits a held that evolves to a spherical wave and at 
distances large with respect to the aperture is effectively a 
plane wave since the phase curvature varies little across the 
aperture. These uncorrelated spherical waves evolve at the 
aperture to angular plane waves with incoherent phases. This 
is expressed in the term e l ^' r where 


25 is the delay due to the arrival of the light from the source. 

In conventional narrow -held Michelson interferometry, 
the action of a Michelson interferometer is to overlay the 
helds from two spatially shifted apertures with an artificially 
induced time delay between one of the helds (delay line) . The 
30 held at the entrance pupil 202 is truncated by the two aper- 
tures, each of diameter D. The center-to-center spacing of 
each aperture is denoted as the baseline B. The relay optics 
208 shifts these two portions of the held such that they are 
overlaid and introduces a time delay via a delay line 206 in 
35 one arm such that they are coherently summed with the inten- 
sity given by: 

/i(*. y,*)= <4 > 


/ 1 ( B x 

B, 

AU 

if B x 

By 

AAi 2 \ 

(14 -T 

,y--,zit- 

" 2\ 

| + £o(*+ — • 

■ y+ T' ; 

b,+ tI) 


where B=(B x ,B. > ,) is the baseline vector, or the spacing of the 
two apertures in the entrance pupil plane and At is the variable 
time delay. The intensity is given by I 1 (x,y,z)=(E 1 (x,y,z;t)l 2 ) 
where 


-> 2 nv 

k = —(or, fa y) 


50 


<■■■) 


i *r 

" 2z;J " 


is the wave vector in the direction of propagation of that plane 
wave component and (a,|3,y) are direction cosines of that 
plane wave component such that a 2 +(3 2 +y 2 =l. For small 
angles, approximations a»sin a, (3«sin (3 and y«l hold. Using 
these approximations 


denotes averaging over the random variable of phase. The 
55 phase is a random variable since each emitter in the incoher- 
ent source emits independently of all other source emitters. 
The intensity is evaluated by inserting the fields from equa- 
tion (3) into (4) and evaluating to show: 


60 


-*■ 2 7rv 2 nv 

k-r = (ax + py ) H z and 


e = e c 


I fax, y, z) = h (x, y, z) + h(x, y, z) + 


(5) 


2R q(Eq(x 



z\t — 


A t\ / B X By 

tFoP -~-,y+-~-,z\ 



yields for equation (1): 


65 


where l l and I 2 represent the intensities over each sub -aper- 
ture. The interference occurs in the cross-term (3 rd term). 
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Inserting the fields (equation 3) into this term and re-airang- 
ing and collecting like terms yields: 

{E 0 E* 0 ) = j j' J\E 0 (a, fa-, v)\ 2 e^ {Bx ? +B y®e^*dadfadvdp ^ 
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where A(x,y) is the aperture function, unity inside and zero 
outside the aperture and (0 X , 0 ) are focal plane positions 
projected on the sky. 

Following the steps outlined above, the interferometer 
equation can be derived. The focal plane intensity, projected 
on the sky, is given by (E 1 +E 2 l 2 )=(E 1 l 2 )+{E 2 l 2 yf 
2 Re (E 1 E 2 * ). Propagating to the focal plane yields: 


Recognizing that I E 0 (a, (3 ;v) 1 2 is the emitted sky intensity at 
each sky angle and frequency provides the interferometer 
equation shown below: 


h(x, y, z) + h (x, y, z) + I 2 (x, y, z) + 


(7) 


2Re{ J J J I 0 (a, fa; v)e~^ l B * a+B y® e ~ 2nvAt dadfadv} 


h(0 X , Oy) = h(9 X , By) 


Hff 


PSF(a -9 x ,fa- 9 y , v)/ 0 (or, fa; v)dadfa 


(E l E‘ 1 ) = J J J PSF(a-9 x , fi- 
B y , v)/ 0 (of, fa\ v)e~ i?J E (aBx+/3B y ) dadfa 


( 11 ) 


Equation (8) below shows the results of integrating all the 
flux in equation (7), i.e. integrate over (x,y): 

h(B x , B y , At) = (8) 

h +/ 2 +2Re{ f f f *>(<*, fa, v)e~ i? T LlBxa+B y^e~ i ^ t dadfadv ) l 


Factoring out l 1 +I 2 =I r yields and noting that both I^BjAt) 
and I 0 ( a ;v) are both real yields the interferometer equation: 


where PSF(a-0 x ,p-0y,v) is the primary beam point spread 
function I 0 (a,p;v)=IE 0 (a,(3;v)l 2 is the sky intensity at each 
wavelength. Equation (11) demonstrates that \ = \ 2 and that 
each takes the form of the sky brightness angularly convolved 
with the primary beam PSF and subsequently integrated over 
frequency and that these two tenns are independent of base- 
line and delay. 

The 3 rd term in equation (11) is the interference term. 
Noting that the time delay can be expressed as a path length 

difference via AL=cAt and by defining: 0 =(0 x ,0 > ,), 

B=(B X ,B > ,) and cx=(a,(3) yields for the focal plane intensity: 


Ii(B , At) = /;r|l + 2 J' J' J " -^-^cos|2tt|-S-c? + vArj|ja(iv| 


(9) 


Equation (9) is one form of the narrow field of view Mich- 
elson interferometer equation. Its significance of is that for 

each baseline position I?=(B X ,B > ,) and each time delay At the 

interferometer observes one real number, ^(3, At), and thus 4Q 
as the baselines are moved and the delay line varied a 3D data 

set in the variables (S,At) is realized. The sky brightness 

function, I 0 ( a ;v), can be estimated at each frequency by 
Fourier techniques. Note that the sky angles have been com- ^ 

bined in a vector a =(a,j5). 

In one aspect, the number of baselines and the maximum 
length of the baselines effectively filter which spatial frequen- 
cies are present in the estimated sky intensity and also limit 
the angular resolution. The length of delay line limits the 50 
spectral resolution and the delay line sampling limits the 
range of frequencies. The field of view is limited to approxi- 
mately +/- 1 X/D where X is the wavelength and D the diameter 
of a single sub-aperture. This narrow field condition on the 
conventional Michelson is why wide field Michelson is 55 
needed. 

Continuing the propagation from the exit pupil 210 to the 
focal plane involves a propagator which steps the field from 
the pupil 210 to focal plane and is given by the standard 
Fraunhofer propagator for each frequency. For an exit pupil 60 
field of the form E P (x,y,v) the focal plane field is given by: 

Ef( 0 x , 0 y , v) = r f A(x, y)Ep(x, y, v)e~ s2R c ( ' 9xX+e y y) dxdy 


/ F (0; B, A L) = 

2 J' PSF(a - 9, v)/o(S; v)jl + cos|27r- (a-B - AL)J jdakf v 


All terms are real, i.e. not complex, and the frequency integral 
extends only over positive frequencies since negative fre- 
quencies do not exist. 

Equation (12) is the interferometer equation which 
expresses the collected focal plane flux as a function of sky 
angle, baseline spacing and path length delay for a given sky 
brightness distribution. This form can be compared with 
equation (9), i.e. the conventional Michelson interferometer, 
and it is seen that the focal plane intensity is now also a 

function of 0 , i.e. the location in the focal plane. The col- 
lected focal plane flux is a function of 5 variables (5 -dimen- 
sional) given by (O^O^B^B^AL). Equation (9) can be 

obtained from equation (12) by integrating over all d 0 with 
the point spread function normalized such that 



PSF(a -0)d0 = 1. 


Equation (12) represents the starting point for estimating a 
hyperspectral data cube. 

Next the disclosure turns to construction of the wide-held 
Michelson interferometer dirty data cube. The interesting 
problem is the inverse problem of constructing the sky bright- 
ness distribution, i.e. I 0 ( a - 0 ), at each frequency or wave- 
length from observations of 1^ 0 ;S,AL). In a wide-held 
Michelson at optical frequencies the delay line is realized by 
moving an optic in one arm of the interferometer in a con- 
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tinuous fashion. The movable section can be in a one or two 
dimensional axis and can be moved linearly or rotationally. 
Metrology is used to determine the relative position of the 
mirror as a function of time and the path length delay is given 
by AL. The baselines are realized by stepping the spacing 5 
between the two collecting apertures in a discrete fashion, i.e. 
the apertures are moved and stopped and then the delay line is 
moved through its range and a focal plane image collected at 
each delay line position, then the baseline mirrors are moved 
and the cycle repeated. This effectively generates the data 

described by equation (12) with a few caveats. Namely 8 is 
not continuous but consists of a finite set of discrete baselines 
and the delay line positions are sampled at discrete times and 
over a finite range and the detector pixels integrate over a 15 
finite area and are also spatially sampling. These effects com- 
bine to give a discretely sampled version of equation (12), in 
addition noise in the detector and instrumental errors also 
causes systematic and random effects in the collection of the 
data and ultimately in the estimation of the image. These 20 
effects are described below but first the basic image construc- 
tion process is derived. 

Let equation (12) be represented as the sum of two inte- 
grals: 


-continued 


^(S-3-AL)]d 


Equation (15) can be converted to its quadrature term by 
evaluating its Hilbert transform with respect to AL to yield: 


H{V(9; B, AL)} = 

— - — I" d v f+ PSF(a - 9 , v) / 0 (5; v)sin[2^- (B-a- AL)1 da 

/..Jflwo J- 00 L c 1 


An analytic form can be constructed from equations (15) 
and (16) as: 


V (a) (0; B, AL) = V{9\ B , AL) - M{v{9\ B, AL)} 


-PT 

'l J 0 J—oo 


PSFld-9, v 


-«2^(£.3-A L) d ^ 


1 f {9 ; B, AL) = 2 f C PSF(a - 0, v)/ 0 (a; v)da dv + '' 1 ' 

do d— CO 

2 £f PSF(a-9, v)/ 0 (3; v)cos[2^(fl-tt - AL)] dadv 30 

Average over delays are denoted as 


I Ave 0) = lim J f ■■■dAL. 

L -+00 LJq 


The average of equation (13) is evaluated as: 


Fourier transforming over AL with v'/c as the conjugate 
Fourier variables yields: 


V ia) {9\ B, v ' ) = — - — PSF(a - 9 , v' )0(a; 

lAve 0 ) J ~°° 


/, (1^) 
v )e c da 


Define the inverse Fourier transform with respect to base- 
line, 8) but with v'/c 0 as the conjugate Fourier variable as: 


V (o) (0, v') = V (a) {9 ; B, v'y^dB 


1 r L _» ^ 

I Ave {9 ) = lim - l P {9\ B, AL)dAL 

Lr+co L do 


= 2 if PSF(a - 9, v)/ 0 (a; v)dadv 
Jo J— 00 

= 2 J^° PSF(9, v) ** /q(0; v)dv 


However only a finite set of discrete baselines 
{8 } _i _ _ N exists, implying that 


n b n+co 

2 r 

j=l-oo J -°° 


B-B i y Z7r T^dB = 


I Ave ( 0 ) is independent of the baseline and is obtained for 

each focal plane position 0 by averaging its interferogram 
over the delay line positions and across all baselines. The 55 
primary beam point spread function is convolved with the 
object and integrated over all frequencies. 

Let the interferograms with the time averaged removed at 
each baseline and focal plane pixel be defined by: 


£ V (a) ( 9 ;Bj; v')e <2 4* 
j= 1 


Inserting equation (20) into (19) and re-arranging gives 
equation (21) shown below: 


v(0; B , AL) = 


I f (9-B, A L)-I Ave {9) 


^Td V r 

i AV e(e) Jo 


- ( + 0(d-v')PSF(d-~9, v')V e-^cM^da 

)d-°° tl 


PSF(d - 9, v)Io(a; v)cos 



11 


US 8,144,331 B2 


12 


Defining 


and recognizing that 


n B /_ 

PSF d (9 , v')£ e~ i 2 ^f 9 5 

as the complex dirty yields: 

10 

V to) (0, v ') = — r°°0(a; v')PSF d (a-9, v')da 
lA Ve 0) J -°° 

2 

= P57^(0, v')**0(9 ; v') 

/a«(3) 15 

where * * denotes two-dimensional convolution. 

Thus the dirty image is given by PSF^( 0 ,v')**0( 0 ,v') 
and thus is given at: 20 


/<*M=^Re{? (o, M} 

25 

This result shows that the dirty image cube, at each fre- 
quency can be constructed from the set of baseline cubes via 
the Hilbert transform, a set of Fourier transforms over each of 
the delay line interferograms, then a discrete inverse Fourier 30 
transform across the set of baselines. 

In the discretely sampled focal plane the sampling is 
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i.e. approximately Nyquist sampled with respect to the diam- 
eter, D, of a single aperture where X=c/v is the wavelength. 
The phasor sum in equation (21) can benefit from sampling of 


A9< 

2\Bj\ 

and since |S / I>D in general, or at least for any useful inter- 
ferometer, direct evaluation of equation (21) may result in 
severe aliasing. One exemplary approach is to employ the 
Fourier shift theorem. Take the Fourier transform of 

y(«)( Q respect to 0 and with iT as the conjugate 

Fourier variable to show: 




is the Fourier transform of V ( ^( 0 jSyv') but shifted to 


v' -► 



in the u domain. Thus, the procedure is to Fourier transform 

V^( 0 ;S,v') at each baseline and frequency into the FT 
domain. Then numerically shift them into a larger grid such 
that the origin is shifted to 



This results in a large grid in the u domain at each frequency 
and each of these is subsequently inverse transformed to yield 

V^(0, v') on a larger sampled grid. This procedure is per- 
formed for each frequency and the results used to estimate the 
dirty image cube. Discrete Fourier transform techniques can 
provide enhanced results with sparsely populated grids. 

l Ave ( 0 ) is sampled at lower resolution than the result from 
equation (23) and hence can be interpolated. This is straight- 
forward since l Ave ( 0 ) does not contain any interference and 
thus its spatial frequencies are limited to the spatial frequen- 
cies of the primary beam. 

The disclosure now turns to a discussion of the process to 
construct a wide-field imaging interferometry testbed (WIIT) 
dirty image cube. FIG. 3 illustrates an example schematic of 
WIIT 300. WIIT 300 is described here from a processing 
point of view. Light from the source module 302, with select- 
able test scenes, reflects off the large collimating parabola 304 
and enters the two movable reflective apertures 306 of WIIT 
300. The source module 302 is at the focus of the collimating 
parabola 304. These apertures 306 move top to bottom on a 
sliding rail and the distance between them is the baseline 308 . 
A series of relay optics 310 fold the beams from each of the 
two apertures around to the beamsplitter 312 which combines 
the two beams in the pupil plane. Subsequently a lens 314 
brings them to focus onto the CCD camera 316. An actuated 
moving carriage which has laser metrology generates an opti- 
cal delay by providing differential path length positions 


n+oo 

—CO 


J-oo j =l 

n b 

= p V (a) (0; Bj, V)e~ i2 ^~ V c B jf do 
j= i 



(24) 
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between the two arms of the interferometer. The source scene 
can be rotated to provide 2D baselines and hence emulation of 
two interferometric dimensional imaging. 

FIG. 4 illustrates an example data format for WIIT data and 
the layout of the WIIT data as it is collected. At a given 5 
baseline position the system 100 collects a data cube. Each 
plane of the cube is a 32x48 pixel image representing the two 
dimensional field of view. In this example, the delay line 
includes about 5000-7000 steps and thus this number of 
frames to each baseline cube. Each baseline cube can be to 
stored as a self-contained Flexible Image Transport System 
(FITS) data cube in Tables format with the metrology data for 
each delay line position. Baseline cubes can be stored in other 
suitable formats as well. In one embodiment, the delay line 
moves in approximate increments of 9 . 89 nm and the detector 1 5 
pixel sizes are 2.475 arcseconds. Other dimensions and dis- 
tances are contemplated. Initial processing loads each base- 
line cube sequentially, parses the differential delay lines posi- 
tions table, and interpolates the data to equal increments of 
the delay such that the delay line sampling is uniform for each 20 
of the baseline cubes. 

FIG. 5 illustrates example WIIT interferograms and vis- 
ibilities. A single frame image 502 of an observed WIIT data 
cube can be a predetermined size, such as 48 pixelsx32 pixels. 
Plots 504A-504F show counts versus delay line position for 25 
the pixel at location 36x16 of the single frame image 502 for 
6 different increasing baselines. The interference pattern 
moves to the right with increasing baseline. An image of the 
visibility amplitudes 506 at the pixel at location 36x16 rep- 
resents the increasing baseline for all 21 baselines (left to 30 
right) and versus delay (top to bottom). The visibilities pat- 
tern moves towards the upper right and diminishes. Visibili- 
ties are obtained via the Hilbert transform. 

FIG. 6 illustrates an example WIIT fractal bar target and 
images used for a test scene. The geometry of the test mask 35 
602 includes 7 bars where light is passed such that each 
successive bar from left to right is closer together by a factor 
of Vi from the previous. The mask, or simulated bar target, 
604 is shown and an observed WIIT image 606, based on raw 
WIIT data, is shown. In one aspect, the left most bars are 40 
resolved while tending to blur together towards the right as 
they pass through the primary beam resolution limit. The 
width of a bar can be 0.85 arcseconds as compared to a 
detector element width of 2.475 arcseconds and is thus unre- 
solved. The simulated dirty image 608 illustrates what the bar 45 
target can look like after processing. 

FIG. 7 illustrates an example image 702 of WIIT one- 
dimensional ultraviolet plane coverage. Imaging of the bar 
target is accomplished by moving the two WIIT apertures to 
21, for example, positions in one dimension resulting in the 50 
UV-plane coverage 704, modulation transfer function (MTF). 
The MTF is completely filled in one dimension except for a 
slight dip 706 near the origin due to the apertures coming 
close together but because of hardware limitations not touch- 
ing. 55 

FIG. 8 illustrates an exemplary WIIT dirty beam for the 21 
baselines. The dirty beam is shown on a linear scale 802, on a 
log scale 804, and as a plot through its center 806. A fringe 
narrow with respect to the core of the primary beam is evident 
in the center and this fringe width defines the limiting reso- 60 
lution of WIIT. 

FIG. 9 illustrates example WIIT spectra for varying base- 
lines. FIG. 9 shows spectra 900 as processed from the delay 
data for the brightest data pixel for each baseline cube. The 
variability shown in the spectra 900 is likely the result of the 65 
processing algorithm and may not perfectly accurately reflect 
actual spectra 900. 
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FIG. 10 illustrates sample WIIT dirty image data results 
and maximum entropy deconvolutions. FIG. 10 shows the 
first estimate of the dirty image 1002 of the one dimensional 
fractal bar target 602 of FIG. 6. The exemplary image 1002 
shown is at one wavelength (600 nm) which corresponds to 
the peak of the spectrum shown in FIG. 9 and the output 
image is 1024x1024 samples. The dirty image 1002 shows 
some edge ringing from the bar targets and the phases of the 
recovered Fourier components of this dirty image 1002 may 
not be entirely correct. The four frames 1004, 1006, 1008, 
1010 below the top frame 1002 in FIG. 10 show a maximum 
entropy deconvolution for 5, 10, 15 and 20 iterations respec- 
tively. The inputs are the estimated dirty image and the cal- 
culated dirty beam shown in FIG. 8 and a specialized form of 
maximum entropy is used which incorporate maximum like- 
lihood constraints and is optimized for ill -posed inversion 
problems such as interferometric aperture imaging. A Hilbert 
transform based algorithm can estimate the visibility ampli- 
tudes and phases in a relative sense across the wide-field of 
view. Maximum entropy methods can assist in better estimat- 
ing the object in the dirty image cube. 

FIG. 11 illustrates an example inverse Fourier transform 
over B with 0 as a conjugate variable. FIG. 11 demonstrates 
the steps 1100 to generate a wide field hyperspectral data cube 
from a set of dirty image cubes. First a system 100, such as 
that shown in FIG. 1, loads the B, cube 1102. The system 100 
removes the mean and renormalizes the loaded cube 1004 as 
shown by the equation inbox 1004. The system 100 generates 
an analytic signal via a Hilbert Transform with respect to AL 
1106. The system uses a Fourier transform to generate inter- 
ferograms 1108 based on the analytic signal. The system 100 
also uses a Hilbert phase transform per B 7 to store the phase 
offsets 1110. The system 100 combines these results to store 
Bvis cubes 1112. The system loops over the set of cubes 
representing discrete sets of baselines 1114 until no cubes 
remain. 

FIG. 12 illustrates an exemplary method embodiment for 
wide field imaging interferometry. The method can be per- 
formed by an appropriately configured system 100 such as 
that shown in FIG. 1. The system 100 first performs the 
following steps for each point in a two dimensional detector 
array over a field of view of an image: gathering a first inter- 
ferogram from a first detector and a second interferogram 
from a second detector (1202), modulating a path-length for 
a signal from an image associated with the first interferogram 
in the first detector (1204), overlaying first data from the 
modulated first detector and second data from the second 
detector (1206), and tracking the modulating at every point in 
a two dimensional detector array comprising the first detector 
and the second detector over a field of view for the image 
(1208). The detectors can be telescopes. The interferograms 
can represent visible light, invisible light (such as infrared), 
sound, radio waves, seismic waves, gravitational waves, or 
any other suitable type of wave. 

The system 100 then generates a widefield data cube based 
on the overlaid first data and second data for each point 
(1210). The system 100 can then generate an image from the 
widefield data cube. In one embodiment, the system selects 
the first and second detectors from a group of three or more 
detectors, such as a grid or array of detectors. FIG. 15, dis- 
cussed below, illustrates on such possible configuration. The 
first and second detectors can be different physical detectors 
for each point. In one aspect, the system 100 gathers more 
than one first and second interferograms at different distances 
between the first and second detectors for each point in the 
two dimensional detector array. In another variation, the sys- 
tem 100 gathers more than one first and second interfero- 
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grams at different locations for the first and second detectors 
for each point in the two dimensional detector array. 

FIG. 13 illustrates an example configuration of two detec- 
tors on a single boom. A single boom 1300, rod, or other linear 
structure attaches to two detectors 1302, 1304. The detectors 
slide back and forth on the boom 1300 to detect at different 
distances from each other. This example configuration pro- 
vides an image of dimensions roughly corresponding to the 
shape and size of the boom 1300. 

FIG. 14 illustrates an example configuration of two detec- 
tors on a rotatable boom 1400. The first detector 1402’ s 
position is fixed on a boom 1400. The boom rotates about the 
first detector 1402 and the second detector slides along the 
boom 1400 at different distances from the first detector 1402. 
For instance, at position 1404, the second detector is at the 
farthest possible distance away from the first detector 1402. 
At position 1406, the boom 1400 has rotated clockwise and 
the second detector has moved slightly closer to the first 
detector 1402. At position 1408, the boom 1400 has rotated 
further clockwise and the second detector has again changed 
distance on the boom relative to the first detector 1402. At 
position 1410, the boom 1400 has rotated almost a complete 
360 ° to position 1404, and the second detector is very close to 
the first detector 1402. With an appropriate number of 
samples from various distances and rotational positions, this 
system can generate a wide-field image equivalent to that of a 
much larger single detector or telescope. 

FIG. 15 illustrates an example grid configuration 1500 of 
more than two detectors. This grid configuration can be static 
or all or part of the grid can be dynamically movable. The first 
and second detector can be different detectors in the grid for 
each iteration or for multiple iterations of the same point. For 
example, the first and second detectors can be detectors 1502 
and 1504. For a next point, the first and second detectors can 
be detectors 1502 and 1506, or 1502 and 1508, or 1506 and 
1510 . In some cases, multiple detectors in the grid can act on 
the same point at the same time. In another aspect, the first 
detector 1502 can pair with multiple second detectors 1504, 
1506, 1508, 1510 simultaneously. The distance between the 
detectors varies through one approach or another. 

Embodiments within the scope of the present disclosure 
may also include computer-readable media for carrying or 
having computer-executable instructions or data structures 
stored thereon. Such computer-readable media can be any 
available media that can be accessed by a general purpose or 
special purpose computer, including the functional design of 
any special purpose processor as discussed above. By way of 
example, and not limitation, such computer-readable media 
can comprise RAM, ROM, EEPROM, CD-ROM or other 
optical disk storage, magnetic disk storage or other magnetic 
storage devices, or any other medium which can be used to 
carry or store desired program code means in the form of 
computer-executable instructions, data structures, or proces- 
sor chip design. When information is transferred or provided 
over a network or another communications connection (either 
hardwired, wireless, or combination thereof) to a computer, 
the computer properly views the connection as a computer- 
readable medium. Thus, any such connection is properly 
termed a computer-readable medium. Combinations of the 
above should also be included within the scope of the com- 
puter-readable media. 

Computer-executable instructions include, for example, 
instructions and data which cause a general purpose com- 
puter, special purpose computer, or special purpose process- 
ing device to perform a certain function or group of functions. 
Computer-executable instructions also include program 
modules that are executed by computers in stand-alone or 
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network environments. Generally, program modules include 
routines, programs, components, data structures, objects, and 
the functions inherent in the design of special-purpose pro- 
cessors, etc. that perform particular tasks or implement par- 
5 ticular abstract data types. Computer-executable instructions, 
associated data structures, and program modules represent 
examples of the program code means for executing steps of 
the methods disclosed herein. The particular sequence of such 
executable instructions or associated data structures repre- 
io sents examples of corresponding acts for implementing the 
functions described in such steps. 

Those of skill in the art will appreciate that other embodi- 
ments of the disclosure may be practiced in network comput- 
ing environments with many types of computer system con- 
15 figurations, including personal computers, hand-held 
devices, multi -processor systems, microprocessor-based or 
programmable consumer electronics, network PCs, mini- 
computers, mainframe computers, and the like. Embodi- 
ments may also be practiced in distributed computing envi- 
20 ronments where tasks are performed by local and remote 
processing devices that are linked (either by hardwired links, 
wireless links, or by a combination thereof) through a com- 
munications network. In a distributed computing environ- 
ment, program modules may be located in both local and 
25 remote memory storage devices. 

The various embodiments described above are provided by 
way of illustration only and should not be construed to limit 
the scope of the disclosure. For example, the principles herein 
can apply to. Those skilled in the art will readily recognize 
30 various modifications and changes that may be made to the 
principles described herein without following the example 
embodiments and applications illustrated and described 
herein, and without departing from the spirit and scope of the 
disclosure. 
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We claim: 

1. A computer- implemented method of wide field imaging 
interferometry, the method comprising: 

for each point in a two dimensional detector array over a 
40 field of view of an image with all steps implemented 
within a computer processor: gathering a first interfero- 
gram from a first detector and a second interferogram 
from a second detector, in a computer modulating a 
path-length for a signal from an image associated with 
45 the first interferogram in the first detector, overlaying 
first data from the modulated first detector and second 
data from the second detector in the computer processor, 
and tracking the modulating at every point in a two 
dimensional detector array comprising the first detector 
50 and the second detector over a field of view for the 
image; and 

generating a widefield data cube based on the overlaid first 
data and second data for each point. 

2. The computer-implemented method of claim 1, the 
55 method further comprising generating an image from the 

widefield data cube. 

3. The computer-implemented method of claim 1, wherein 
at the first and second detectors are selected from a group of 
three or more detectors. 

60 4 . The computer- implemented method of claim 1 , wherein 

at least one of the first and second detectors is different for 
each point. 

5. The computer-implemented method of claim 1, the 
method further comprising gathering more than one first and 
65 second interferograms at different distances between the first 
and second detectors for each point in the two dimensional 
detector array. 
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6. The computer-implemented method of claim 1, the 
method further comprising gathering more than one first and 
second interferograms at different locations for the first and 
second detectors for each point in the two dimensional detec- 
tor array. 

7 . The method of claim 1 , wherein the gathered first and 
second interferograms represent visible light. 

8. The method of claim 1, wherein the gathered first and 
second interferograms represent a wave other than visible 
light. 

9 . The method of claim 1 , wherein the detectors are tele- 
scopes. 

10 . A system for wide field imaging interferometry, the 
system comprising: 

a processor; 

a module configured to control the processor to perform the 
following steps for each point in a two dimensional 
detector array over a field of view of an image: gathering 
a first interferogram from a first detector and a second 
interferogram from a second detector, modulating a 
path-length for a signal from an image associated with 
the first interferogram in the first detector, overlaying 
first data from the modulated first detector and second 
data from the second detector, and tracking the modu- 
lating at every point in a two dimensional detector array 
comprising the first detector and the second detector 
over a field of view for the image; 

and a module configured to control the processor to gen- 
erate a widefield data cube based on the overlaid first 
data and second data for each point. 

11 . The system of claim 10 , the system further comprising 
a module configured to control the processor to generate an 
image from the widefield data cube. 

12 . The system of claim 10 , wherein at the first and second 
detectors are selected from a group of three or more detectors. 

13 . The system of claim 10 , wherein at least one of the first 
and second detectors is different for each point. 
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14 . The system of claim 10 , the system further comprising 
a module configured to control the processor to gather more 
than one first and second interferograms at different distances 
between the first and second detectors for each point in the 

5 two dimensional detector array. 

15 . The system of claim 10 , the system further comprising 
a module configured to gather more than one first and second 
interferograms at different locations for the first and second 
detectors for each point in the two dimensional detector array. 

to 1 6 . A tangible computer-readable storage medium exclud- 

ing all transitory storage media storing a computer program 
having instructions to control a processor for wide field imag- 
ing interferometry, the instructions comprising: for each 
point in a two dimensional detector array over a field of view 
15 of an image: 

gathering a first interferogram from a first detector and a 
second interferogram from a second detector, modulat- 
ing a path-length for a signal from an image associated 
with the first interferogram in the first detector, overlay- 
20 ing first data from the modulated first detector and sec- 
ond data from the second detector, and tracking the 
modulating at every point in a two dimensional detector 
array comprising the first detector and the second detec- 
tor over a field of view for the image; and 
25 generating a widefield data cube based on the overlaid first 

data and second data for each point. 

17 . The tangible computer-readable storage medium of 
claim 16, the instructions further comprising generating an 
image from the widefield data cube. 

30 18 . The tangible computer-readable storage medium of 

claim 16, wherein at the first and second detectors are selected 
from a group of three or more detectors. 

19 . The tangible computer-readable storage medium of 
claim 16 , wherein at least one of the first and second detectors 
35 is different for each point. 



